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Abstract 

In order to investigate the features of the classical approximation at high temperatures 
for real time correlation functions, the plasmon frequencies and damping rates were 
recently computed numerically in the SU(2)+Higgs model and in the pure SU(2) the- 
ory We compare the lattice results with leading order hard thermal loop resummed 
perturbation theory. In the broken phase of the SU(2)+Higgs model, we show that 
the lattice results can be reproduced and that the lattices used are too coarse to ob- 
serve some important plasmon effects. In the symmetric phase, the main qualitative 
features of the lattice results can also be understood. In the pure SU(2) theory, on the 
other hand, there are discrepancies which might point to larger Landau and plasmon 
damping effects than indicated by perturbation theory. 
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The dynamics of non-abelian gauge fields at finite temperature has attracted some 



theory. However, some physical quantities of interest are non-perturbative and the only 
known method of calculating them on the lattice is the classical approximation [![]. 

One problem of the classical approximation is related to the fact that the high 
momentum modes do not decouple from the dynamics of the low momentum modes. 
In the naive perturbative expansion, momenta of order T lead to corrections in the 
Green's functions which are proportional to T 2 . For Green's functions with soft external 
momenta |p| ~ g 2 T, where g is the gauge coupling, these "corrections" can be as large 
as or even larger than the tree level Green's functions. These large corrections have 
been named hard thermal loops (HTL) and a consistent perturbative expansion requires 
that they be resummed 0. 

In the resummed perturbation theory the gauge field propagators reflect two distinct 
collective phenomena. One is the plasma oscillations which involve simultaneous os- 
cillations of the low and the high momentum degrees of freedom in the plasma. The 
characteristic frequency of these oscillations is proportional to gT. Secondly, there is 
an effect called Landau damping which is related to an energy transfer from the low 
momentum degrees of freedom to the high momentum ones. 

In the classical approximation the physics of the high momentum degrees of free- 
dom is not correctly described. The hard thermal loops correspond to classically UV- 
divergent contributions ||. How the hard thermal loops affect non-perturbative corre- 
lation functions of the low momentum fields is an open problem. 

The physical picture sketched above is based on 1-loop resummed perturbation the- 
ory. It is therefore interesting to see whether it also shows up in non-perturbative 
lattice studies. In a recent paper |ll]] the correlation functions of two gauge invariant 
operators in the SU(2)+Higgs model (in continuum notation), 



were computed with classical lattice simulations. Here (p denotes the Higgs doublet, 
r a is a Pauli matrix, $ is the matrix $ = (iptp) with ip = ir 2 (p*, and V denotes the 
space volume. In the broken phase the correlator of W?(t) is given by the gauge field 
propagator. It was claimed that the plasmon frequency corresponding to this correlator 
is independent of the lattice spacing a. This contradicts the picture above since in the 
classical theory the plasmon frequency is proportional to T/a instead of T 2 . We will 
demonstrate here that the qualitative features of the correlators can be understood 



attention recently |l ] [16]. A lot of work on this subject has been done in perturbation 
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using perturbation theory and that the claim of ref. |TT| concerning the broken phase 
plasmon frequency is not justified. 



Ambj0rn and Krasnitz |L2[ have reported a measurement of the gauge field correlator 
in the Coulomb gauge in the symmetric phase. We will see that many of the features 
observed can be qualitatively understood in perturbation theory, but quantitative dis- 
crepancies remain. 

Real time correlators and the HTL effective action 

The plasmon properties of the classical Hamiltonian SU(2)+Higgs gauge theory could 
be computed by solving the classical equations of motion perturbatively with given 
initial conditions, and by then averaging over the initial conditions with the Boltzmann 
weight |[|. However, it appears that in perturbation theory it is simpler to perform 
the computation in the full quantum theory and to take the classical limit h — > only 
in the end. We use this approach. 

In the classical limit the operator ordering in a correlation function is irrelevant. 
Thus the quantum expression of which we are going to take the classical limit can be 
anything but a commutator. This still leaves several possibilities, and we choose to 
consider 

C Q (t, P) = \J d 3 xe'^ x (0(t, x)0(0, 0) + O(0, 0)O(t, x)> . (3) 

In a perturbative computation, one first computes the corresponding two-point Green's 
function Po(ia> n ,p) in Euclidean space for the Matsubara frequencies u n . Then one 
performs an analytic continuation which depends on the real time correlator in question. 
The symmetric combination in eq. (0) is given by 



—e- lwt [1 + 2n(u)} Im V a (iu n ^u + te,p), (4) 

-oo 2lX 



where 



<^ = e/3 J_ x ( 5 ) 

is the Bose distribution function. 

To evaluate the w-integral in eq. (^) numerically, it is convenient to rewrite it (for 

t > 0) as an integral in the upper half of the complex tu-plane. In this way one 
avoids integrating along the poles and discontinuities on the real cj-axis. Denoting 
uj = to\ + iu 2 , one gets for t > in the classical limit n(u) — ► T / (Huj), 

/OO (jr.)-. 

— 7 — e- i ^ t \imV (coi + iw 2 ,p). (6) 
-oo 2m Id i + lid2 
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Note that this expression is independent of the imaginary part co>2 > 0. The second 
term on the r.h.s. of eq. (|]) vanishes for the equal time case t = 0. A further virtue 
of eq. @ is that one can take the limit p —>■ inside the integrand since u>2 > 
guarantees that there are no singularities on the integration contour. In the remainder 
of this letter we will discuss only classical correlation functions^ and we will therefore 
omit the superscript 'classical' from the l.h.s. of eq. 

We thus have to evaluate Vo(w,p). As discussed above, a consistent perturbative 
expansion requires the resummation of the hard thermal loops In other words, 
the high momentum modes p>T (or p > a -1 in the classical theory) are integrated 
out, giving an effective theory for the low momentum modes p<gT. This amounts to 
including the hard thermal loops in the tree level action. Denoting P = (u n ,p), the 
quadratic part of the Euclidean HTL effective theory in momentum space is 



r( 2 ) _ 1 \a aci 
^HTL — 2 V " 



P 2 5 „v _ pl x p v + n^ L (p) 



P 2 ~ \ m2 H + S HTL 



(7) 



where mu is the Higgs mass. Eq. (0) gives directly the gauge and scalar field free 
propagators T>°^(uj, p), V{uj, p), to be used in the computation of Voiyj, p). Note that 
E H tl is momentum independent. 

Of the HTL self-energies, we will need especially the components I1^ TL (P), S H tl- 
The expressions are || H, |10| , after the replacement u n — > —iu, 



Ktl(P) = -9 2 (2N + N s )hJ- 



d 3 q dn{uj q ) v l v^u 



(27r) 3 d,LU q LU — p l V l 



where, for the SU(2)+Higgs model, N — 2, N g — 1. Furthermore, 

du), 



-<i 



(10) 



dqi 

where u q is the tree level dispersion relation; for a space-time with discretized spatial 
dimensions (a cubic lattice), 



-v&w^), (ii 



auj q 



Eqs. (|S|), (|S|) hold actually for a generic dispersion relation |ID| . In ref. ||, scalar 
electrodynamics was considered?] for which case g 2 (2N + N s ) has to be replaced with 
2e 2 in eq. (§). 

3 Except in eq. ( |37| ) and in the discussion thereafter. 
4 The result for n^ TL given in [Bj has the wrong sign. 
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Note that in the theory of eq. (0), the gauge field A is included. In the classical lat- 
tice simulations fTI, I^fl, in contrast, one puts A = and imposes the Gauss constraint 



explicitly. 

The hard thermal loop self-energy simplifies in the limit p = 0, which is relevant for 
the operators in eqs. (|l|), (Q). In that limit, the gauge field self-energy is 

n°° TL = o, n°v L = o, ng TL = ^u4, (12) 

where u)w is the plasmon frequency, 

^^' 2 (1 a) 



to 



Let us also introduce the "scalar plasmon frequency" 

1 

= ~2 m ^ + Ehtl ' ^ 
In the continuum limit of the full quantum theory, for which |v| 2 = 1, one has 

= g^2N + N,)~, (15) 

<& = 4 m « + ( 6A + 3 j2 )l5S- (16) 

whereas in the classical limit n(uj) — > T / (fiu) on a cubic lattice, one gets 
Here 



2 n \ 4" J 4tt a 



E = 1 r /2 d s x = 3.1759H535625 (19) 

7T 2 J-tt/2 X)t s i n x i 

is a constant which can be expressed in terms of the complete elliptic integral of the 
first kind . 



The basic problem of the classical real time simulations can now be expressed as 
follows. For the scalar field, one has a mass parameter in the classical SU(2)+Higgs 
theory which does not break gauge invariance. It can be tuned such that the high 
momentum modes produce the correct quantum expression in eq. and that the 

lattice spacing dependence in eq. ( |TE| ) is canceled. For the gauge fields, in contrast, the 
local classical SU(2)+Higgs theory does not allow a mass term and hence the divergent 
classical expression in eq. flT7| ) cannot be arranged to coincide with the quantum ex- 
pression in eq. fli"5|). This divergence should also show up in the classical simulations. 
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This problem is specific for time dependent correlation functions: in the static case 
only the ui n = sector requires a resummation and IIjj TL (0, p) = 0. 

Finally, let us fix some notation. The continuum field ip in eqs. ([[]), (|2|) is related to 
the dimensionless lattice field tp in [|TTJ by 



<p*<p 



4 



G 



T(3g t ™ ag2T 
The operators measured on the lattice are then 

1 



H 

W? 



W 3 X 

1 



'N 3 



^^ f (x)^(x), 

x 

^'Tr [$ t (x)C/ i (x)$(x + e i a)r a 



(20) 

(21) 
(22) 



where £7j(x) is the link operator. The correlators measured with these operators are 
denoted by Cg(t), ^S ab SijCy^(t), and they are related in the continuum limit a — > 
to the correlators C#(i), C^At) of the continuum operators in eqs. ([]]), (||) through 
C H (t) = a(3lT 2 C R (t), C$ tij (t) = ^dija-^^Cwit). The factor 3 in the definition 
of Cw{t) corresponds to a sum over the different isospin components. 

The broken phase 

Let us first consider the broken phase. The real time correlators of the operators H(t), 
W^(t) have been determined in the classical approximation in the broken phase of the 



SU(2)+Higgs model by Tang and Smit JTT 



We parameterize the scalar field in the broken phase as 



$=^b(T) + 0o+ir°0 o ] 



(23) 



Then, in continuum notation, the lowest order terms in the gauge invariant operators 
of eqs. ([]]), (@) become 



H(t) 

wt(t) 



const. + 
1 



d 3 xv(T)(f) (ti x), 



d s xv(T) m w (T)A$(t,x) - di<i> a (t,x) 



(24) 
(25) 



where mw{T) = gv(T)/2. Due to the integration Jd 3 x, these operators correspond to 
zero spatial momentum, p = 0, so that the second term in eq. (^5|) does not contribute. 
It follows that the leading terms in the correlation functions Cnif) and C^^(t) are 
given by the propagators of the elementary fields 0o and A^. 
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To compute the required propagators, one can set p = in the time dependent part 
of eq. (|6]), so that eqs. fll3|), (|T4"D can be used. At zero spatial momentum the full 
Euclidean (Matsubara) propagators for the fields A", O are of the form 

5 ab 5- ■ 1 

where *II, *£ denote the parts of the self-energies which are generated radiatively 
within the HTL effective theory. The tree-level terms are 

{u h w f = m 2 w (T) + 0.05379/3 G (/T) 2 , (27) 
(^) 2 = m%(T), (28) 

where mw{T), m#(T) are the mass parameters generated by the Higgs mechanism in 
the static theory: mw(T) = gv(T)/2, rrijjiT) = u% + 3\v 2 (T). Here the scalar mass 
parameter was tuned to its correct value by using the counterterm of the static classical 
theory and the fact that Ehtl is momentum independent. For (^{y) 2 , in contrast, 
the /^-dependent part comes from eq. (|17|) and cannot be removed. Parametrically, 



~ g 2 T 2 and TI(u;]^), *S(u;^) ~ g 3 T 2 , so that the dominant contributions to the 
H and W" 1 plasmon frequencies (appearing as C(t,p) ~ C exp(— Tt) cos(a;pit + 5) in 
the correlator) are just u\ and according to eq. (H). The damping rates Y are 
related to *n(c^) and *S(a;^). 

In Fig. m we compare the leading order plasmon frequencies in eqs. (^3), (^8|) with 
the lattice results of ref. [[11]]. The value of v(T) has been determined from the 1-loop 
effective potential. The zero-temperature parameters used in |TTJ correspond to run 
= 80 GeV. It is seen that the lattice results are remarkably close to the leading order 
perturbative results. In particular, we conclude that the gauge field plasmon frequency 
uj\y diverges in the continuum limit according to eq. (|27|) , while the scalar plasmon 
frequency uj\ remains finite and equals the static screening mass at leading order. One 
sees that the lattice is so coarse that it is difficult to notice the divergence of since 
this is shadowed by the finite mw(T). Thus one is in a sense not close enough to the 
continuum limit. It should also be noted that the amplitude of C^At) dies out as 
\/{l)w) 2 in the continuum limit. 

The damping rates of the gauge and scalar fields are parametrically of order g 2 T 
and therefore, in contrast to the plasmon frequencies, they are classical. This has 



been demonstrated explicitly in a scalar field theory [T^, 14]. A full computation of 



the damping rates in the broken phase of the SU(2)+Higgs theory is missing at the 



moment. However, the order of magnitude can apparently be understood []TT| using 



the known symmetric phase gauge and Higgs elementary field damping rates |18|, [19] 
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Figure 1: A comparison of the leading order perturbative plasmon frequencies with 
those determined on the lattice in [IT]]. The continuous curves are the perturbative 
values from eqs. (|27|), fl28|), and the squares are the data points from Table 5 in |TT] 



Open symbols correspond to ui^y and filled to u^. The corresponding {3 G values are 
given next to the symbols. 



The symmetric phase 

Let us now turn to the symmetric phase. The correlators of the composite operators in 
eqs. (0), (0) have been determined in the symmetric phase of the SU(2)+Higgs model 



by Tang and Smit [11]. In addition, the gauge field correlator of the pure SU(2) model 
has been measured in the Coulomb gauge by Ambj0rn and Krasnitz fL2"|. 



Consider first the composite operator correlators measured in pT| . In the symmetric 
phase, the composite operator character of H(t) and W^(t) manifests itself more clearly 
than in the broken phase. For instance, the leading term in is 



1 



/ N 3 



•o(x + eia)(f) a (x.) - o (x)0 a (x + e^a) - e afec (/> fe (x)0 c (x + e^a) , (29) 



which does not contain the gauge field A" at all. 

Due to the fact that no gauge fields are involved, the leading terms in Cnit) and 
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K + P 

Figure 2: The lowest order diagram contributing to Cuif) and C^^{t) in the symmetric 
phase. The solid lines denote scalar propagators. 

Cyyjj(t) can be easily computed: both are given by diagrams of the type depicted 
in Fig. |2|. To evaluate them one has to start with a Matsubara external momentum 
p = iu n . Then the sum over the loop frequencies is written in terms of an integral in 
the complex k plane (see, e.g., [pOj] ). Only then can one continue to arbitrary complex 
values of Pq. For the operator 0(t, x) = (p*<p(t, x), the diagram in Fig. |^ finally gives 



ImP^Po + ie, P 
d 4 k 



4 



(2tt)< 



n(k ) - n(k + p )\ lmV(k + ie, k) lmV(k + po + ie, k + p). (30) 



This expression^ is equivalent to a corresponding one derived in ref. [^l]]. The correlator 
C H (t) = l(H(t)H(0) + H(0)H(t)) for the operator H(t) of eq. <g) is given by C H {t) = 
C ' 9 i v {t, p = 0), where CLt„,(t, p) is obtained from eq. (|J). 

To get the leading contribution in Cu{t), one uses the free propagators for which 



lmV {vcc (k + ie, k) 



7T 



2 aA 2 + ^ 



5 (*„ - yW^) - ^ (^0 + 



Then the integrals over k and po can De performed in eqs. (|3 
p — > 0, one obtains 

d 3 k 1 



C H (t)=T 2 J 



1 + cos 2Jk 2 +uj 2 H t 



(2tt) 3 (k 2 + u 2 H ) 2 
On the lattice (in the continuum limit) this corresponds to 



(31) 



In the limit 



(32) 



C s (t) 



dx- 



x 



4vr 2 ^7o ~~(x 2 + z 2 ) 2 
where z = aun/^ and i = t/a. 



1 + cosiAtVx 2 



(33) 



5 Eq. ( pp[ ) can be written in other forms by a change of integration variables, but then one may get 
a wrong result for the free case, if e is put to zero inside the integral and eq. (pi]) is used. 
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Figure 3: The leading order perturbative correlator Cg(t) in lattice units in the 



symmetric phase at T/T c — 1.52. To be compared with Fig. 8 in The oscillation 



frequency is 2luh (this agrees at leading order with the static screening mass). 

The vector correlator Cyfjjif) has an additional factor of kikj in the integrand com- 
pared with eq. Q3"i2|). Therefore the continuum limit of C^^it = 0) does not exist. On 
the lattice one obtains 



Cw(t) = C d3x iZ^7 Xt 2 + Z ^ t 1 + C ° S ( 4 ^sin 2 x, + ^) 



(34) 



The correlation function Cfj(t) is shown in Fig. |3| and C^(t) in Fig. ^. These should 
be compared with Figs. 8, 10 in flU , respectively. It is seen that the qualitative features 
can be understood quite well with the leading order results. Note, in particular, that 
the scalar correlation function C^{t) in Fig. ([$]) is oscillating with an amplitude which 
decreases with time. It should be emphasized that this decrease is not related to 
damping (remember that Fig. [3] shows the tree level results without any interactions 
and damping occurs only through interactions). The decrease is rather due to the 
fact that there is a continuous spectrum of frequencies lo > 2u H causing a destructive 
interference in the phase space integral, eq. (§^). This shows that it is difficult to 
determine a damping rate from the gauge invariant operator Cg(t). 
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Figure 4: The leading order perturbative correlator C\y(t) in lattice units in the 
symmetric phase at T/T c = 1.52. To be compared with Fig. 10 in |TTJ. The oscillation 
period here is proportional to the lattice spacing a. 



Let us then try to estimate how higher order corrections could modify the qualitative 
behavior of C#(t). Consider the effect of self-energy insertions in the scalar propagators 
T>(ko, k). The self-energy has an imaginary part. Therefore the scalar propagator does 
not have poles at ko 

k 2 + u> 2 H , 
One can 



±y k 2 + uj% (the lowest order result for Cn{t) is due to these 
poles). Since the imaginary part of the self-energy is small compared with 
Cn{t) will nevertheless still be dominated by the region k w ± 
therefore approximate the scalar propagator as 

1 



k 2 + 



H- 



-(k + iT k ) 



V(k + ie,k) 



where the width T k is given by (see, e.g., |pl| 
Inserting this into eq. (|30~1) and using eq. 



+ k 2 



2 • 
H 



(35) 



T k = -ImE( 
, we find 



k 2 + ^,k)/(2 v /k 2 + ^ 



C H (t)=T 2 j 



d 3 k 



-2F k t 



(2 7 r)3[k 2 + ^][k 2 + ^ + r 2 



1 + cos 2 v /k 2 



10 



ii 



t - 2a h 



(36) 
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where = arctan yTk/yk 2 + uj 2 H y Thus the effect of damping should show up in 
Fig. |]such that the constant part (corresponding to the first term in the square brackets 
in eq. (0)) decays away at large times. This is indeed the qualitative behavior observed 
in m. 



Finally, we consider the transverse gauge field correlator CV(£,p) measured in [12 
Gauge fields can only be defined in a particular gauge, which in |L2] was chosen to 
be the Coulomb gauge. We let the external momentum point into the ^-direction, 
p = e 3 (2nk/L), where L is the spatial extent of the lattice and k is an integer. Then 
CA t (t, p) is given, e.g., by the correlator of Af, 

C At (t,p) = ±y>ze^"(A?(t,xK(0,O) +A?(0,OK(i,x)>. (37) 

Let us first recall some features of C A ,(t, p) and of the corresponding analytic Green's 
function J) A t (k>, p) for |p| ~ g 2 T in the quantum theory. After the HTL resummation, 
V A .(uj,p) has poles at uj ~ ±oj w , where u>w is given by eq. (f[5|). These poles lead 
to an oscillation of C At (t, p) on the time scale \fh/{gT). In addition, "D^ t (o;,p) has a 
discontinuity on the real tu-axis for uj < |p|. This discontinuity is related to Landau 
damping and it gives a contribution / Lan dau(£, p) to C At (t, p). The function / La ndau(£, p) 
does not involve any oscillations and just constitutes a decaying background for the 
superimposed plasmon oscillations. The time scale on which /Landau (^p) varies is 

Higher order corrections to V At (uj, p), which are not included in HTL effective action 
but are generated radiatively within that theory, lead to a damping of the plasmon 
oscillations. The plasmon damping rate T is of order g 2 T and has been computed for 
p = in ref . |18| . 



These two different damping effects manifest themselves in the correlator C At (t, p) 
in quite different ways, so that its functional form is expected to be 

CU t (*> P) ~ ^exp(-rt) cos(u w t + 5) + landau (^, p)- (38) 

The time dependence of /Landau (^ p) becomes non-perturbative for t>l/{hg A T) 0, 
where /Landau (^ p) is expected to vanish. This means that CA t (t, p) can be computed 
perturbatively up to a non-perturbative constant C At (0,p) as long as t 1/(% 4 T). 

In the classical lattice gauge theory one expects a similar qualitative behavior. In the 
order of magnitude estimates for the quantum theory one has to replace h — > Ta. The 
analytic structure of V At (u,p) is more complicated than in the quantum case. The 
HTL resummed V At (oo, p) depends not only on the magnitude but also on the direction 
of p pi. In particular, there are directions of p for which there is a discontinuity for 
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arbitrarily large values of uo ||10|| . Nevertheless, the qualitative behavior of CA t (t, p) 
should be given by eq. 



Unfortunately, the numerical results for CA t {t,p) in |TjJ have been normalized to 



(7^(0, p) which cannot be computed in perturbation theory Comparing our per- 
turbative estimates with the non-perturbative results therefore requires some model 
assumptions about C*A t (0,p). 

In order to account for the plasmon damping effects, we include the leading order 



damping rate T in the HTL resummed propagator V^(u,p). Since the momentum 



dependence of T is not known we use its value for p = which is T = 0.176g 2 T |T 
We can then write for the transverse components in eq. (|37D, analogously to eq. (|35|), 

5 ab 

^^.p) = _ (u)+<r) , +p , +nBrLKp) - (») 

Consider first the zero external momentum case, p = 0. Eq. (^) only makes sense 
when r 2 < p 2 + Ren^rpL^, p). In the latter term in eq. (@) there should be no 
problem, since for p = 0, Il^ TL (u;,p) — > u^y, but for the first term (i.e., for the static 
limit t = 0) the inequality is not satisfied (remember that Hh TL (0, P) = 0). One might 
try to regulate the real part of the gauge fixed self-energy phenomenologically with 
a "magnetic mass"; letting M = yrra 2 ^ + T 2 + u; 2 ^, M = Jrr^^ n + ujyy, it follows 
from eq. (||) that this would give 

C^O) = m2 T +r2 . (40) 

magn ' 



Cjl it) ~ Cju (0) = -— - + — ^e~ Tt cosf Mt - arctan ^J) . (41) 
tV ; tV ; M 2 MM V MJ 

Parameterized this way, one can indeed find quite reasonable agreement with Fig. 6 



in 



12[ , but only if T is chosen to have a large value, T ~ (0.8 . . . 1.0)g 2 T. Otherwise one 
is getting too many oscillations in CV (t), not observed in |T^]. For these large values 
of T, the magnetic mass parameter is in fact favored to be small or zero. However, 
these fits are quite phenomenological, and thus we will not consider them any more. 

For p ^ 0, the t — -part of the leading order correlator is still parametrically 
non-perturbative, but at least it is formally finite for r — > so that the perturbative 
approximation might be numerically reasonable. To get a feeling about the momentum 
scales in question, note that for the value $q = 14 considered in |^2| , |p| = 0.69g 2 T x k 
and u>w = 0.87g 2 T (the latter can be obtained from eq. (l27f) with m w {T) = 0). We 
have evaluated numerically both the HTL self-energy in eq. (§) and the remaining uj\- 
integral in eq. for k = 1, . . . , 4. The resulting correlators are shown in Fig. |5|. This 



should be compared with Fig. 6 in [12 
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Figure 5: The gauge field correlator in the pure SU(2) theory. To be compared with 



Fig. 6 in |12[]. Here t is expressed in lattice units: t = t/a. Unless otherwise stated, 
Pg = 14. 



The main effects to be seen in Fig. |^ are the following. For k = 1 one can see 
an oscillation in Fig. j|, in contrast to Fig. 6 in |12| . Thus one could say that the 
non-perturbative plasmon damping rate is larger than the perturbative estimate T = 
0.176<7 2 T. Indeed, one has to go to a much larger value, T ~ 1.0g 2 T (> |p|, u)\y), to get 
enough plasmon damping. Another observation to be made at k — 1 is that in the 
correlator is already very close to zero at t/ (5q = 0.10. This is not quite so in Fig. |, but 
one has to wait much longer for the correlator to vanish. Thus the non-perturbative 
Landau damping effects seem also to be larger than the leading order perturbative HTL 
result. At k = 2 there is no large difference in Landau damping any more, but the 
plasmon damping appears to be somewhat too weak even with T ~ 1.0g 2 T. Finally, 
at k = 4 one gets reasonable agreement between the perturbative and lattice results, 
provided that T ~ 1.0g 2 T. One can also see that the Pg~ dependence is reproduced; 
the plasmon frequency thus diverges according to eq. (JT7|). 

We have also tried the expression — uj 2 — 2iuT in the denominator of eq. (|39|) , so 
that the real part is not modified by Y 2 . The qualitative conclusions and the preferred 
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value T ~ 1.0g 2 T remain the same. Based on those curves, one would nevertheless say 
that even at k = 2 there is too little Landau damping in the perturbative estimates, 
but at k = 4 one again gets good agreement. 



Summary and Conclusions 

We have computed several quantities related to real time correlation functions in the 
classical SU(2) and SU(2)+Higgs models on the lattice, using hard thermal loop re- 
summed perturbation theory. 

Our results for the gauge field and scalar plasmon frequencies in the broken phase 
are in remarkable agreement with the numerical lattice simulations in ref. [[U]]. We 



have reiterated that the classical gauge field plasmon frequency is divergent in the 
continuum limit and we have demonstrated that this is consistent with the results of 



ref. [II], where the plasmon frequency was claimed to be lattice spacing independent. 

For the symmetric phase we have computed gauge invariant scalar and vector corre- 
lators as functions of time at the lowest order in perturbation theory. Furthermore, we 
have estimated the effect of higher order corrections. Our results are in good agreement 



with ref. JTT|]. We have shown that it is difficult to extract damping rates from the 
measurement of these correlators. 

Finally, we have studied the correlator of the transverse gauge field in pure SU(2) 
gauge theory. While the qualitative features of the numerical simulations in ref. |12] are 
consistent with our perturbative estimates, there appear to be significant quantitative 



discrepancies. The damping of the plasmon oscillations observed in [12| appears to be 
much stronger than one would expect from the perturbative result for the damping 
rate [18|. This is puzzling because the damping rate is of the order g 2 T and should 



therefore have a classical continuum limit. However, one should keep in mind that the 
perturbative estimates are reliable only if the plasmon frequency is much larger than 
the damping rate which is not the case for the lattice spacing used in ref. [[12] . 



Acknowledgments. We are grateful to E. Berger, P. Overmann, O. Philipsen, 
M.G. Schmidt and I.O. Stamatescu for useful discussions. 

Note added 

After this paper was submitted, we were informed by J.Smit that in the revised ver- 
sion of Ref. [[ll|] , Tang and Smit have weakened their conclusions concerning the lattice 
spacing independence of the W plasmon frequency in the broken phase, so that their 
statements are now in better accordance with ours. We thank J.Smit for communica- 
tion on this issue. 



In view of the interpretation in 22 1, let us stress that the discussion around eq. (|35|) is 
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not meant to be a consistent quantitative estimate of the higher order corrections. We 
just want to see in which qualitative way the higher order corrections might manifest 
themselves. 
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